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c3 Abstract 



The evaluation of BPS Wilson loops in M = 6,D = 3 Chern-Simons matter theory is reduced to 

ordinary matrix integrals via localization technique. It is easy to check that the vacuum expectation 

' _(' value of 1/2 BPS Wilson loops at leading order in planar limit agrees with the regularized classical 

JL string action, via AdS/CFT. Then the subleading terms in principle can be calculated by treating 

r~| the string theory semi-classically In this article we calculate the one-loop determinant for the 

^ fluctuation modes of holographic Wilson loop as IIA string m the dual geometry AdS, x Cp3. The 

^L fermionic normal mode frequencies are expressed in terms of the hypergeometric function, and we 

3^ compute the one-loop effective action numerically. The discrepancy with localization formula is 

due to the zero mode normalization constant, which is yet to be determined. 
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I. INTRODUCTION 

Wilson loops are essential objects in the study of gauge field theories. In the context of 
the AdS/CFT correspondence pQ, they have dual description as a macroscopic fundamen- 
tal string [21 E]. In this article, we are mainly interested in the M2-brane conformal field 
theory as Chern-Simons matter model, suggested by Aharony, Bergman, Jafferis and Mal- 
dacena (ABJM) [4J. The supersymmetric Wilson loop operators in ABJM model, with dual 
geometry AdS^ x CP^, are studied earlier in [SjiSj. The computation of their expectation 
values can be greatly simplified if one utilizes the localization technique [H]: when we put 
the gauge theory on S^, the full path integral is reduced to an ordinary matrix integral [TU] . 
It is a fascinating achievement that at strong coupling the free energy scales as N^/"^ and the 
coefficient is related to the internal space S'^ , precisely as predicted by AdS/CFT ^Tj, [T2] . 

According to the matrix model calculation at strong coupling and planar limit, the 1/2- 
BPS circular Wilson loop's vacuum expectation value is (up to a framing-dependent phase) 

where A is the 't Hooft coupling constant. On the other hand, the gravity side computation 

from classical string solution is e '^. The next-order correction for S* = — In(iy) should be 

ln2 ;^ 0.69, and it is our goal to see if this number can be reproduced as one-loop correction 

on string world sheet. 

From the fluctuation lagrangian around 1/2-BPS holographic circular Wilson loop, we 

find that the string one-loop determinant is given as 

^_r ^ det(-V^-l)det^(-V| + ^) 

det(-V2 + 2)det^(-V2) ' ^ ' 

It turns out that part of the fermionic normal mode frequencies in the numerator are given 

in terms of hypergeometric functions. This is in contrast with the Wilson loop of IIB string 

in AdS^ X 5^, where the frequencies are logarithm of rational functions and the sum is given 

exactly using the Gamma function |il3j. We evaluate F numerically and extract the finite 

piece after regularization, and obtain F^eg ~ —1.1. 

In Section [TT] we setup the notation and calculate the quadratic lagrangian for string 



fiuctuation around 1/2-BPS circular Wilson loop. In Section III, we calculate the normal 



modes and discuss how their sum can be regularized numerically. In Section IV we discuss 



how to resolve the discrepancy between field theory and supergravity side results. 



II. OPEN STRINGS AND THEIR FLUCTUATION LAGRANGIAN 

We consider type IIA open strings in AdS^ x CP^ background which preserve | super- 
symmetry. This geometry is conjectured to be dual to A/" = 6, Z^ = 3 Chern-Simons field 
theory [1] with U{N) x U{N) gauge symmetry and levels {k,—k). In the convention we 
adopt here, the D = 10 supergravity solution takes the following form. 

F2 = fcJcps, F4=^VoU,54- (3) 

o 

Rs sets the length scale of this background, the metric tensor ds\^g^, ds'^^-j, are scaled to 
have radius one, and J^ps represents the Kahler 2-form of the internal space. The AdS/CFT 
correspondence relates the string and Chern-Simons description in the following way. 



RsHa' = {2T,'\Y'\ (4) 

where A = N/k is the 't Hooft coupling constant. For simplicity we will henceforth set 

k = l. 

It is convenient for us to use Poincare coordinates for AdS space, 

X 
ds\dSi = -^(-dt"^ + dr^ + r^dcp"^ + dz^) . (5) 

For a circular Wilson loop with radius 1, a simple solution is given as 

z = Vl — T^ . (6) 

In conformal gauge r = 1/ cosh a and the induced metric on worldsheet is 

dsl^ = ^{da'^ + dr"^), < a < oo, < r < 27r . (7) 

sinh a 

Note that this is hyperbolic with scalar curvature R^"^^ = —2. In order to regularize the 
divergence of the classical action, we introduce a cutoff at 2; = e or equivalently at a = eo 
which are related via e = tanheo. The regularized value of the classical action is [SHT] 



So = -Rl/a' = -V2^. (8) 

Now we are to consider the fluctuation modes around this classical solution. Similar 
computations have been performed in a number of articles including [T3HT9] . For the bosonic 



sector, the computations should be very similar to those of the circular Wilson loop in 
AdS^ X S^ presented in [13]. One easily finds that after the gauge fixing, there are two 
modes from AdS space with effective mass parameter 2, and there are six massless modes 
from CP . Altogether they account for the denominator of pi). 

For the fermionic part, up to quadratic order the ^-symmetric Green-Schwarz action is 

Sp = -^ [ dM^h-'S'' - e'-'S'^WpaDi^'e'' . (9) 

zna J 

where 5^"^ = diag(l,— 1), pa = T^daX^^ E^^. X^^ parametrizes the ten-dimensional space- 
time, Ta is gamma matrix, E^j is vielbein, and hat is the worldsheet metric. The spinors 9^ 
and 6"^ have opposite chirality, i.e. 

Tii^i = e^ , Tiie^ = -e^ . (lo) 

The covariant derivative for spinor field is spelt out as [H] 

D'f = (^da + \daX''utfTAB'j ' ^^.X^^i^i/^BCr^^las)^^ (11) 

After some calculation one can rewrite the fermion fiuctuation lagrangian simply as 

C = i^lC^, IC = Vh{T'Vi-tTy4Toi). (12) 

We note that this expression is obtained after rotating the spinor by a unitary matrix 

S = exp I — Fis I , tan a = - . (13) 

\1/ also satisfies P+\E' = \1/ with P+ = (1 + roirii)/2, d = 2 gamma matrices r* satisfy 
{ri,Tj} = 2hij, and 

13/4 = ^ (3123 + (r45 + Tgy + r89)rnroi) . (14) 

It is obvious that r3/4 is hermitian and traceless. When diagonalized, it can be written as 
for instance diag(l, 1, 1, 0) ® diag(l, 1, —1, —1). This implies that we should have 4 massless 
fermionic modes, and 12 modes with mass 1, on the worldsheet. We note here that this 
result is in agreement with similar analysis done for instance in [TTj [T9] . 



For the computation of the determinant, we might as well consider the square of Dirac 
operator. We consider 

R(2) 

A^ = {^T'V^ + Ts/,Toif = -V| + ^ + r^/4 , (15) 

where V^^ = -^Vi^y/gg^^Vj) and for the solution we have here R^'^^=—2. 

Our results so far can be summarized in the following expression for one-loop partition 
function for fluctuation modes. 

^_ det^/^(-V|-i)det^/^(-V| + ^) 

det2/2(-V2 + 2)det6/2(-V2) ' ^ ^ 

Note that in the denominator V^ is the usual scalar Laplacian, while V^^ is understood 
to contain spin connection for spinor fields. One can repeat the same computation for a 
straight line which is also 1/2-BPS and we have checked the result is again given exactly as 



(16). 



III. CALCULATION OF THE DETERMINANT 



Now let us consider the evaluation of (16). Thanks to the axial symmetry of the string 



worldsheet, we can easily perform the mode expansion for r variable. We impose periodic 
(anti-periodic) boundary condition for bosonic (fermionic) fields. Then Z can be expressed 
using determinants of ordinary second-order differential operators. More concretely, we have 
for instance 

det(-V2) = Y[det[smh'' cr{-dl + n^)] (17) 

nez 

det^(-V|) = Yl det[sinh^(T(-9^ + z/2 + |cothV + z/cothcr)] 

X det[sinhV(-9^ + z/2 + lcothV-z/coth(T)]. (18) 

The conformal factor sinh a cancel between bosonic and fermionic determinants. We define 

<i = In [det(-a2 +n^ + 2cschV)/C] , (19) 



B3 



U. 



UJ^,' 



uP 



\n[det{-dl + n^)/C], (20) 

In [det{-dl + z/^ + z/ cothcx + \ coth^ a + |cschV)/C] , (21) 

In [det(-9^ + iy^ + iy coth a + \ coth^ a - |cschV)/C] , (22) 



where we have included C = det(— (9^) as an overall normalization. The 1-loop effective 
action can be written as 



i y: (^r+c^i+scf +3.^^) 



r^-lnZ = 5^(a;f + 3a;f)-- ^ (c.- + c.- + Scf + Sc.^). (23) 

n£Z iyeZ+1/2 

It turns out that each sum ^ Un is divergent and there is an ordering problem. This 
problem is of course commonplace in quantum field theory, and for the energy correction 
of spinning strings in 74^5*4 x CP'^ the ordering issue has been addressed in [T71 HH 120] . 
Here we follow the prescription in [131 EI] : one introduces a regulator fi in the process of 
synchronizing the summation indices for bosonic and fermionic modes. For small /i, we have 



r.e. = $:^"""'(^n' + 3^n') 



neZ 



1 J2 e-^l^l^-^^ 
uez+i/2 



<'+.."+ 3c." + 3^5) 



neZ 



,-Fl 



,-Fl 



.Fl 



e-M|n| ( 4^m ^ ^2u;B3 _ ^f]_^^^ _ ^^^^^^ _ ^^._^^^ _ ^^^^^^^ 



q, ,F3 q, ,F3 q, ,F3 

"'J'^n+1/2 ~ <J^n-l/2 ^ 'J^-n-1/2 






+ (e-MM _ e-'^l«+V2|)(^^^^^ + ^^i_^ + ^^Fs + 3^^3_^ 



1/2 



_l_ (e"^l"l _e"''l""^/^l'ii'''''^^- 



1+1/2 



F3 



I/2J 



.F3 



)(^n-l/2 + '^-n+1/2 + 3^^-1/2 + 3W-n+l/2) 



n=0 

Here we have defined 

1 

2 
1 

2 



Go 

Gn 



{2u^' + Qu^"" - uj(l - uj^_\,^ - 3u;f/| - Sa;^?^) 



4a;f ^ + 12wf ^ - u^L ,o - wf_li lo - uj—_, i^ - uj 



^n+l/2 



■'n-1/2 



■^-n.-l/2 



Fl 

-n+1/25 



q, ,F3 q, ,F3 q, ,F3 q, ,F3 

-OCJ„+l/2 — <J^n.-l/2 ~ <J^-n-l/2 ~ 'J^-n+1/2 



(n>0) 



G' = lim^y. 



'fm 



n>0 






,i^l 



,F1 



■'n+1/2 



n-1/2 '^n-1/2 ^^-71+1/2 



+3w„_^l/2 + 3W-n-l/2 ^ 3a;n-l/2 ~ ^^-n+1/2 



(24) 



(25) 



(26) 



(27) 



A. Calculation of the frequencies 



To evaluate u^ and co^, following p^l22] - [2l] we utilize the Gelfand-Yaglom theorem: For 
a differential operator O with periodic boundary condition in a G [a,b], the product of all 



eigenvalues can be alternatively obtained by solving the homogeneous differential equation 
Oip = with initial condition tp{a) = ipo{a) = 0, ip'la) = ipQ{a) = 1. In particular, 

detC -^{b) 



detCn 



MbV 



(28) 



where Oq = ~d^. For our problem originally a ranges in < cr < oo, but we will introduce 



both UV and IR regulators and consider instead tQ < a < L. Eq.(28) will be used for 



non-zero modes, while for the zero-modes we take Neumann boundary conditions at L, and 



we need to use ip'{b)/ipQ{b) instead on the right hand side of (28). 

For the bosonic part the operators are exactly the same as the counterpart in AdS^ x S^ of 
IIB string theory, and we simply import the results in [131 . For large L (eo is not necessarily 
small yet.). 



eM^n) 



exp^o;. 



B3\ 



(|n|+cotheo) |n|(L-£o) 
2|n|(|n|+l) e 

cotheo, n = 

g|n|(L-eQ) 

1, n = 



n^O 



nj^O 



(29) 



(30) 



Let us now turn to the fermionic modes. For the differential operators associated with 
fermionic fluctuations, we find it useful to introduce a new variable 



( = coth cr. 



(31) 



Note that for < a < oo, we have 1 < C < oo. We start with the equation associated with 
oj^^. We should originally consider 



{-dl + iy'^ + iy coth a + \ coth^ a + |cschV)?/^(cr) = 0. 
The two linearly independent solutions can be chosen as follows (for u ^ ±^] 



(32) 

(33) 
(34) 



Writing down the solution with appropriate initial condition and taking the limit L — )■ cxd, 
we obtain the following result. 



w„ 



.Fl 



In 


g{>.+l/2)(L-eo) /i+^ 
21/+1 Y 2e2 


, ^ 


In 
In 






/2 

/2. 


4i/2 


-1 £2 


V ^+A 



(35) 



z/< -1/2 



Here we introduced e = tanheo for cutoff of 2;-coordinate. z/ = ±| are studied separately, 
and in particular z/ = — 1/2 is the fermionic zero mode and we have used Neumann boundary 
condition. 

For the other fermionic determinant oj^^, we may employ the following reparametrization 



Ha) = {c + ir^'-'^\c-ir''^'-'^'y{c). 



Again the differential equation is easily solved, and we choose the basis 



u,(a) = (c + ly/'-'/'ic - l)-/2-i/4 

v,{a) = (C + ly/'-'/'iC - l)-/2-V4 



Co (^ + ^y^-'^' 



dx. 



(36) 

(37) 
(38) 



Except for z/ = —1/2 which is zero-mode, the frequency is then (before taking L — )■ oo limit) 



a;f = In 
And in the limit L — )■ oo, 

fcotheo / _ j^Nz/-l/2 



/•cotheo / _ 1 ^1^-1/2 



(39) 



dx 



5(e-2^o;z/+ 1/2,0) z/ > 1/2, 
-^e"'"'^''''^'^ z/<-l/2. 



(40) 



and one should substitute this into (39). Here we have expressed the integral in terms of 
the incomplete beta function, 

B{x;a,b)= f t^'^^l-tf-^dt. (41) 

Jo 

Since z/ is half-integer for our purposes, we may do the integration explicitly and obtain ^ 






2Fi(n, l;n + l;x). 



(42) 



For V = —1/2 we study separately with Neumann boundary condition. Summarizing, we 
have 



c.:^ = < 



In 
In 
In 



e(^+l /2)(L-eo) / 2, 



l/2)CL-£o) 2e 17/1 111 

,+1/2) J—,-2Fi{n,l-n + l-j 



l^e^ 



2e 

-(^+l/2)(Z,-Eo) / 2. 



-2(i/+l/2) V 1+e 



u > 1/2 
z/ = -1/2 
z/< -1/2 



(43) 



-'^ It is also a special case of the Lerch ^-transcendent, i.e. 2Fi(n, 1; n + 1; x) 



n$(x, l,n). $ is defined as 
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B. The regularized action 



We are now ready to go back to (24) and evaluate the finite part. First of all, one can 
easily convince oneself that 

G" = 2(L-eo), (44) 

so G' should have no contribution after regularization. 

It turns out that (?„ — ?■ for large ra, but the series ^Gn for given e is logarithmically 
divergent. We do the sum for large A and drop terms proportional to In A. After rather 
tedious but straightforward computation, we may rewrite ^ G„ in the following way. 



Eg.. 



n=0 



bri 



T 

J- n. 



2 
In 

In 



In 



-4L 



2^(1 + 6)e 

:i- 6)2(1 + 36) 



-(A-l)ln 



1 + e 



+ JE^"-|E^"' (45) 



n=2 



n=0 



{n + \)\n + l){n - I) 



2e 



/n(e)/„+i(e) 



(46) 
(47) 



Here we have introduced a shorthand notation /„(e) = 2-^1(^5 1; ra + 1; ^3^)- 'S'„, T„ are chosen 
such that they converge to zero as n — )■ 00. 

One can see that the total sum is independent of cutoff L, as it should be the case. It 



is also obvious that the finite part of the first term in (45) is ^ln2. The large-A behavior 



of ^ Sn can be studied using Stirling's formula. After we drop the terms proportional to 
In A, l/e,lne etc. 



m 



ln(327r). 



(48) 



reg 



For the summation of T„, unfortunately we are not able to find the sum in closed form 
for large A. We will resort to numerical methods. Our strategy is as follows. We first fix e 
and consider A — ;■ 00. Since the sum is log-divergent, 

A 

^ T„ = /(e) In A + g{e) + (subleading in A). (49) 

n=l 

We can read off f{e),g{e) from a least-square fit after evaluating the sum numerically for a 
number of large values for A. To obtain a regularized value, we now concentrate on g{e). 
This function is also divergent as e — t- 00, and we find it is very closely approximated by 



9 



(6) 



FIG. 1: Dots represent ^(e) from least-square fit of numerical sum ^ T„ against /(e) In A + g{e). 
We used values A = 2500, 2510, •• • ,3500. Solid line corresponds to g{e) = 0.8182 - 1.039^ - 
l.Oll'^. We calculated for 61 points in 0.07 < e < 0.13 and find x^ = 1-68 x 10"^. 

g{e) = a + (3— + J- + 5 log - in the leading orders. The numerical results versus this curve 
is shown in Figure [TJ When we implement this method however, one has to be careful since 
the result depends rather sensitively on the choice of cutoffs A and e. It is not surprising 
since the series is not convergent, after all. We want to send e — )■ eventually, but since 
we take A — ;■ oo first, e should not be too small, i.e. eA ^ 1 should be always satisfied. 
We have tried different ranges for A, e until the result for ^ ^ S^ is reasonably close to the 
analytic result. For 2500 < A < 3500 and 0.07 < e < 0.13, our numerical result is —2.2432, 
where the exact value is — | ln327r = —2.3052. We use the same values of A, e to evaluate 
^T„, and obtain the final result 

Treg = -1.106 (50) 

IV. DISCUSSION 

Treg should bc Compared to the field theory result ln2 = 0.6931, and certainly the dif- 
ference is not negligible. Recall that for 1/2-BPS Wilson loop in AdS^ x S^, there was 
also a discrepancy and it was deemed to come from the normalization of the zero modes 
[T3t [25| [26] . As far as we know, this coefficient is not determined for AdS^ x CP^, let alone 
AdS^ X 5*^ in Type IIB. We note that the normalization convention of holographic Wilson 
loops in ABJM theory was discussed in Section 5.2 of 



10 



To bypass the normalization problem and check the validity of string one-loop compu- 
tations, we may study other supersymmetric Wilson loop operators and calculate the ratio 
between physically different BPS Wilson loops. In the field theory description, there are 
1/6-BPS Wilson loops with {W) ^ ^/\J2 e^^{V2^) . While 1/2-BPS Wilson loops are 
pointlike in CP^ and break the global symmetry SU{A) into ^[/(S), 1/6-BPS ones preserve 
only SU{2) and it is natural to expect that they are smeared over CP^ G CP^ Oil]. We plan 
to construct such classical string solutions in AdS^ x CP^ explicitly and study its fluctuations 
in a separate publication. 



Although we only studied circular Wilson loops in detail here, (16) is valid for a straight 
line as well and our computation is easily extendable to the D = 3 analog of quark- ant iquark 
potential calculation using holography. Of course in principle a similar computation in 
string theory side can be checked against the localization calculation for any supersymmetric 
Wilson loops. Let us emphasize that recently a two-parameter family of string solutions 
interpolating the circle and a pair of straight line Wilson loops in A^ = 4, D = 4 super 
Yang-Mills theory was studied in [271 EH] • It is also pointed out that the angle dependence 
of general BPS Wilson loop operators can be related to interesting physical quantities such 
as cusp anomalous dimension, radiation emitted by a moving quark etc. [29H33] . For a 
recent study of cusp anomalous dimension in AB JM model, see [19] . With such applications 
in mind, it will be intriguing to construct general BPS Wilson loops and pursue their exact 
evaluation. 
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